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Abstract 

Global distribution of aeolian dust is simulated from 1981 to 1996 with the Goddard Ozone 
Chemistry Aerosol Radiation and Transport (GOCART) model. The results are assessed 
with in-situ measurements and the Total Ozone Mapping Spectrometer (TOMS) aerosol 
products. The annual budget over the different continents and oceans are analyzed. It is 
found that there is a maximum of 25% difference of global anuual emission from the min- 
imum in 1996 to the maximum in 1988. There is a downward trend of dust emission over 
Africa and East Asia, of 6 and 2 Tg yr— 1, respectively. The inter-annual variability of dust 
distribution is analyzed over the North Atlantic and Africa. It is found that in winter most 
of the North Atlantic and Africa dust loading is correlated with the North Atlantic Oscil- 
lation. The GOCART model indicates that a controlling factor of such correlation can be 
attributed to dust emission from the Sahel. The Bodele depression is the major dust source 
in winter and its inter-annual variability is highly correlated with the NAO. However, it 
is not possible to conclude without further analysis that the North Atlantic Oscillation is 
forcing the inter-annual variability of dust emission and in-turn dust concentration over the 
North Atlantic. 
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1 Introduction 


Mineral dust has an impact on the different components of the Earth system. Sev- 
eral studies have shown that dust particles, by absorbing and scattering solar radia- 
tion, modify the atmospheric radiative budget (e.g. Tegen et al. (1996); Sokolikand 
Toon (1996); Weaver et al. (2002)). Also, it could play a positive role in reducing 
global warming by greenhouse gas C0 2 . The carbon fixation by phytoplanktons in 
the oceans acts as a sink for C0 2 . Aeolian dust deposition is the primary source 
of bio-available iron in the iron limited open oceans and effectively control phyto- 
plankton blooming ( Martin and Gordon, 1988). Another important effect of dust 
particles is their role in the photochemical production of tropospheric ozone by re- 
ducing by as much as 50% the photolysis rates (e.g. Dickerson et al. (1997); Liao 
et al. (1999); Martin et al. (2002)) and by providing reaction sites for ozone and ni- 
trogen molecules (e.g. Prospero et al. (1995); Dentener et al. (1996)). Finally, dust 
particles affect air quality ( Prospero , 1 999) and are potential vectors for long range 
transport of bacteria ( Griffin et al., 2001). It is thus important to better understand 
the long term variability of dust distributiona and what can be the controling factors 
of such variability. 

In-situ measurements from the 60’s have shown strong daily, seasonal and inter- 
annual variations of dust concentration over the Atlantic ( Prospero , 1 999). Aerosol 
satellite data have been helpful to locate the major dust sources ( Herman et al. 
(1997); Prospero et al. (2002)) and to study the variability of aerosols distribution 
with data since the 70’s ( Herman et al. (1997); Torres et al. (2002)). Also dust 
sources are located on all continents, North Africa seems to be the most productive 
with most dust plumes transported to the North Atlantic. In this paper, we will focus 
on dust emission fron North Africa and dust distribution over the North Atlantic, 
although they will be compare with emissions from other continents and deposition 
over other oceans. Hurrell (1995) have shown that the circulation and precipitation 
over Europe and the North Atlantic is modulated by the North Atlantic Oscillation 
(NAO) with a period of about 8 years. Moulin et al. (1997) found a weak correlation 
between the NAO and the long term variability of in-situ of dust concentration as 
well as with satellite data. Transport models solve explicitly the emission, transport 
and removal processes. Therefore, they can provide more detail information on the 
major processes controlling the spatial and temporal variability of dust distribution. 
Ginoux et al. (2001) have developed a transport model driven by assimilated me- 
teorology which can reproduce successfully dust seasonal variation at the global 
scale. This paper is focusing on the year-to-year variation of the budget over dif- 
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ferent continents and oceans, and on its correlation with the NAO. After describing 
briefly the transport model used to simulate dust distribution, the budget over the 
continents and oceans will be discussed before comparing the model results with 
in-situ and remote sensing data. Finally, the effects of the North Atlantic Oscillation 
on dust emission, concentration and mass column will be analyzed. 


2 Model Description 


The model uses for this study is the Goddard Ozone Chemistry Aerosol Radiation 
and Transport (GOCART) model. The GOCART model simulates the distribution 
of dust, sulfate, carbonaceous (organic and black carbon), and sea-salt aerosols. 
Each model component has been described in details elsewhere ( Chin et al. (2000); 
Ginoux et al. (2001); Chin et al. (2002)). Here we will briefly describe the dust 
component of the GOCART model. 

Dust size distribution is calculated by solving the continuity equation for a dis- 
crete number of size bins. The size distribution of mineral dust extended from 0. 1 
/jm for clay to several hundreds /z m for sand. However, the volume (or mass) of 
particles larger than 10 /zm is several orders of magnitude lower than for smaller 
particles, and the maximum of the volume is around 2 /zm radius ( Dubovik et al.,, 
2002). The size distribution for particles larger than 1 /zm is primarly constrained by 
gravitational settling. In GOCART model, the size distribution is discretized into 
four size bins: 0.1-1, 1-2, 2-3, 3-6 /zm radius. The physical parameters of these 4 
bins are given in Table 1 . The continuity equation includes macroscopic advection 
by winds, parameterized eddy diffusion and moist convection. The removal mech- 
anisms include dry deposition at the surface by impaction, wet deposition in and 
below clouds, and gravitational settling. Dust is uplifted by wind over preferential 
sources which have been associated with topographic lows. The model has a hori- 
zontal resolution of 2° latitude by 2.5° longitude and 20 vertical sigma layers from 
the surface to 1 mb (~ 50 km). All processes are driven by assimilated meteoro- 
logical fields by the Goddard Earth Observing System Data Assimilation System 
(GEOS DAS) at NASA Goddard ( Schubert et al., 1993). 

2. 1 Dust Sources 


The approach used in GOCART to identify the major dust sources is based on the 
analysis by Prospero et al. (2002). Using the Total Ozone Mapping Spectrometer 
(TOMS) aerosol index (http://toms.gsfc.nasa.gov), they have identified and char- 
acterized the geo-morphological nature of the major dust sources. Based on this 
analysis and the previous work by Herman et al. (1997), a global dust-source func- 
tion has been defined as the probability of sediments accumulated in topographic 
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depressions with bare surface ( Ginoux et al., 2001). Figure 1 shows the global dis- 
tribution of the source function on one dgree grid, and the climatological TOMS 
absorbing aerosol index. The correspondence of the maxima of the source function 
and the TOMS index, in the arid and semi-arid regions, is striking. 


2.2 Dust Emission 


Dust uplifting into the atmosphere is mainly initiated by saltation bombardment 
(sand blasting). In GOCART, the vertical flux of dust particles is assumed to be 
proportional to the horizontal flux of sand particles, and it is approximated by an 
expression similar to the empirical formula developed by Gillette and Passi (1988): ’ 


F P = 


C S SpU^Q m (U\Q m Ufp) if UlOm > Uf 

0 otherwise 


( 1 ) 


where C is a dimensional factor equal to 1 pg s 2 m -5 , S is the source function 
described by Ginoux et al. (2001), uio m is the horizontal wind speed at 10 m, u tp is 
the threshold velocity for class p , and s p is the fraction of each size classes given in 
Table 1. 


The threshold velocity for wind erosion is calculated from Iversen and White ( 1 982) 
numerical formulation with the simplications proposed by Marticorena and Berga- 
metti (1995). The expression of the threshold velocity for class p, in units of m s -1 , 
is given by 


u tp = 0.13 



6x 10~ 7 
PpgQp 5 


^1.7638(4.6 x 10 6 [$i- 56 + l] 0 092 



( 2 ) 


where p p is the particle density (kg m~ 3 ), g is the gravity (9.81 m s -2 ), is the 
effective diameter of the class p (m), p a is the air density. The values of p v and 4> p 
are given for each 4 classes p in Table 1 . To take into account the bonding effect of 
soil moisture, the expression 2 is modified as in Ginoux et al. (2001), 


W twp 


u tp x (1.2 + 0.2 log 10 w) if w < 0.5 
oo otherwise 


( 3 ) 


where w is the soil moisture which varies from near zero values in arid region to 1 
for water. 
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3 Model Results 


The dust distribution is siumlated from January 1981 to December 1996. The in- 
stantaneous dust concentration for each 4 bins at every model grid point is archived 
every 6 hours. 


3. 1 Global Budget 


Figure 2 shows the global annual budget for 16 years. It includes dust emission and 
removal by wet and dry deposition. Although the lifetime of dust particles is about 
1 week, the annual budget is not equilibrated: the annual deposition by dry and wet 
deposition is systematically higher than the emission. This is an inherent problem 
of transport models which cannot conserved mass because of time interpolation. 
The error is of the order of 1%. Dry deposition contributes for 90% of total dust 
removal. The ratio between dry and wet removal varies by less than 10% over the 
16 years simulation. The maximum difference of annual emission is about 23% 
with the lowest emission in 1996 (1950 Tg) and the highest in 1988 (2400 Tg). 

Table 2 gives the budget over 5 different regions: North Africa, South Africa, North 
America, South America and Asia. The annual emission from North Africa is 
around 1400 Tg yr -1 which represents 65% of the global emission, while Asia 
contributes for 25%. The contribution from the other regions is relatively low and 
of the order of the standard deviation of dust emission from North Africa. Over 
the 16 years simulation, there is a downward trend of dust emission from Africa 
and Asia of 6 Tg yr -1 , and 2 Tg yr -1 . The deposition rates follow the same trends. 
The difference between the annual rates of emission and deposition gives a mea- 
sure of dust export from the region. The only region with a negative difference, 
import from other regions, is North America with about 30 Tg yr -1 which is 3 
times the amount emitted from this continent. The ratio between dust export and 
emission gives a measure of the efficiency of source regions to provide dust nutri- 
ents the ocean biosphere. The highest efficiency ratio is for Australia (0.3) followed 
by North Africa (0.2), and the lowest is for Asia (0.1), excluding North America. 
The major factors affecting the efficiency ratio is the proximity of the sources to the 
ocean and the altitude and speed of the dust plumes. In Australia, the major dust 
source is in the Lake Eyre region which is about 600 to 800 km from the Indian 
Ocean. On the other hand, the most active dust source in Asia (in the Taklamakan 
desert) is located at more than 3000 km from the Pacific Ocean. Table 3 gives the 
annual dry and wet deposition rates for 6 oceanic regions. The highest deposition 
rates are in the North Indian Ocean while the values in the South Indian Ocean is 
a factor ten lower. The highest contribution of wet deposition is over the North Pa- 
cific (50%) and the lowest is over the South Atlantic (10%). Such large difference 
can be explained by the travel time from the dust source to the ocean. The dust par- 
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tides deposited in the South Atlantic are coming from Patagonia (Chile) which is 
bordered by the Atlantic Ocean and from the Makgadikgadi pan (Botswana) which 
is about 1000 km from the sea shore. 

Table 4 gives the winter, summer and annual budgets over Sahara (North of 2 1 ,25°N) 
and Sahel (South of 2 1 ,25°N) in North Africa. The emissions from Sahara are twice 
the corresponding emissions from Sahel, but ther standard deviations are equiva- 
lent. This means that there is a stronger inter-annual variability of dust emission 
over the Sahel. In Sahara, the emissions are 30% higher in summer than in winter, 
while in Sahel the emissions are slightly higher in winter. In Sahara the wet deposi- 
tion is negligible compare to the dry deposition, but in Sahel it varies strongly with 
season: in winter it is negligible but in summer it represents 22% of the removal 
rate. This can be explained by the movement of the Inter-Tropical Convergence 
Zone (ITCZ) which occupies its northernmost position in summer. South of the 
ITCZ, the monsoon carries precipitating clouds. From Table 4, it appears that the 
inter-annual variability of annual and seasonal emission or deposition rates is dif- 
ferent between Sahel and Sahara. The maximum rates in Sahara date for all seasons 
to 1984, while in Sahel the corresponding year of maximum varies between sea- 
sons: 1983 for winter, and 1988 for summer and annual mean. For both regions, the 
year 1996 has the lowest annual emission. For both regions, there are equivalent 
downward trends of dust emission. 


4 Comparison with In-situ and Remote Sensing Data 

4.1 Monthly Climatology of Surface Concentration 


Figure 3 shows the comparison between the observed and simulated monthly sur- 
face concentration at 16 sites. The results are very similar to the results presented 
by Ginoux et al. (2001), although they have used a parameterization of u t which 
was not appropriate for fine dust but for sand particles. As before, the model is per- 
forming correctly in the dusty regions but overestimate dust concentrations in the 
remote regions of the North hemisphere: over the Atlantic at Mace Head and over 
the Pacific at Midway. In the South hemisphere, the model results are within the 
standard deviation of the measurements. 


4.2 Inter-annual Variability of Concentration 


The surface concentration has been measured for more than 20 years at Barba- 
dos and Miami by the University of Miami ( Prospero , 1999), and for a decade at 
Izana by the New Mexico State University ( Arimoto et al., 1995). Figure 4 shows 
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the comparison of the measured and calculated surface concentrations at Izana, 
Barbados and Miami. The correlations between simulated and observed concentra- 
tions are 0.69 for Izana, and 0.79 for Barbados and Miami. The correlations for the 
monthly anomalies are lower but still significative (99.9% probability of correla- 
tion), they are 0.47, 0.56, and 0.3 1 for Izana, Barbados and Miami. 


4.3 TOMS Aerosol Index 


The TOMS instrument on board the satellite Nimbus 7 measured aerosol backscat- 
tering radiances at 340 and 380 nm from November 1978 until May 1993. In this 
study, we use the data from the period 1981-1990. These measurements cover the 
Earth daily with a 50 km resolution at satellite footprint. The local passing time is 
about 1 1:30 am. By taking the difference between the measured and calculated ra- 
diances for a purely molecular atmosphere, Herman et al. ( 1 997) defined the TOMS 
Aerosol Index (TOMS AI) as follow 

-*»&)}• (4) 

where I m is the backscattered radiance measured by TOMS at the given wave- 
length and I c is the radiance calculated using a radiative transfer model for a pure 
Rayleigh atmosphere. The TOMS AI is a qualitative indicator of the presence of UV 
absorbing aerosols. An inversion procedure that retrieves aerosol properties from 
the TOMS radiances has been developed by Torres et al. (1998). They have shown 
that for absorbing aerosols, like dust particles, the optical thickness can be derived 
from the TOMS AI knowing the values of single scattering albedo, the altitude of 
the aerosol layer, and the surface pressure. Here, an aerosol index is calculated from 
the satellite viewing angles and the simulated size distribution. The methodology 
consists to first calculate the optical thickness at 380 nm from the relation: 

r 380 = ^2 Tk = Y] 7 — ■ — (=1 for k < 4 and k = i for i > 4 (5) 


where r* is the optical thickness at 380 nm for 7 bins k, Qk is the extinction ef- 
ficiency at 380 nm, A/j is the column mass loading of bin i of the 4 transported 
size bins, a* is the fraction of each 7 sub-bins, rk is the effective radius, and p pi is 
the mass density of the size class i. The values of Q*(380nm ) are calculated using 
Mie theory and using the real and imaginary parts of the refractive index derived 
by Colarco et al. (2002). Second, the single scattering albedo at 380 nm ( 0 ^ 330 ) is 
calculated by the relation 


Ylk=l ^k't'k 

W 380 — 

7380 


(6) 
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where cu* is the single scattering albedo of particle radius r * and is calculated using 
Mie therory and the same values of the refractive index as Qk . The values of r^, a*, 
Q/t. and Lu fc are given in Table 5, and the values of p pi are given in Table 1. Third, 
the centroid of mass Z is calculated by the formula 


Z = 


iCijAzj 

Ei Mi 


(7) 


where zj is the altitude above ground at level j among the nlev model levels, Cij is 
the mass concentration of size class i at level j, Azj is the thickness of the grid cell 
at level j, and Mi is the mass column. With r 38 0 , w 38 o, Z, and the TOMS instrument 
viewing angles, an aerosol index is calculated using a look-up table. This table 
has been built using radiative calculations ( Torres et al., 1998) for a limited set of 
values of r380, a; 380 , Z, and viewing angles. The intermediate values are linearly 
interpolated. 

Figure 5 shows the comparison between the observed and calculated aerosol index 
during a dust storm from Sahara in March 1988. Around the 27 March 1988, a 
dust plume was produced from Sahara and transported by the Azores high towards 
Europe. Two days later, the plume is separating in two branches when West and the 
other East. The West branch formed two days later an half circle over most of the 
North Atlantic. The model reproduces the complex pattern of the plume as well as 
the amplitude of the aerosol index. 

The aerosol index has been calculated for the 10 years simulations from January 
1981 to December 1989. Figure 6 shows the global distribution of the correlation 
coefficient between the observed and simulated aerosol index. Over the arid regions 
the correlation coefficient is higher than 0.8, as well as over the North Atlantic. 
Over the North Pacific there is a 15 degrees band with significant correlation but 
most Asian dust plumes are moving eastward at higher latitudes. It is important to 
understand that the TOMS aerosol index includes all absorbing aerosols including 
dust and black carbon, while the simulated aerosol index considers only dust par- 
ticles. So, if the seasonal and inter-annual cycle of the simulated dust distribution 
is in phase with other absorbing aerosols there will be a positive correlation. Such 
phase correlations are apparent in the burning regions in South America. When 
comparing surface dust concentrations at Midway, we showed that the simulated 
values are overestimated. On the other hand, vertical profiles indicate that the peak 
of dust concentration is in the upper troposphere, and not at the surface. The TOMS 
aerosol index increases almost linearly with the plume altitude ( Hsu et al., 1999). 
So, even with relatively small amount, the absorption of solar radiation by aerosols 
in the upper troposphere can produced higher index than larger amount near the 
surface. Additional information would be needed to confirm the presence of dust in 
the upper troposphere over the subtropical Pacific. 
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5 Effects of North Atlantic Oscillation on Dust Distribution 


In this section we examine the inter-annual variability of dust distribution in rela- 
tion with the North Atlantic Oscillation (NAO). The NAO exerts a strong influence 
on the large-scale variations of both the atmospheric circulation and the hydrolog- 
ical cycle in the Northern hemisphere. Hurrell (1995) have defined a NAO Index 
which is calculated by taking the difference between normalized sea-level atmo- 
spheric pressures between Lisbon, Portugal, and Stykkisholmur, Iceland. Winters 
with high NAO indices are characterized by a deepening of the Icelandic low as- 
sociated with a stronger Azores anticyclone. This yields to higher surface pressure 
and drier conditions over Northern Africa. During low NAO conditions, there is an 
increase of precipitation over the Mediterranean and North Africa. Moulin et al. 
(1997) have shown that both pattern and intensity of the transport of African dust 
are affected by the NAO. More recently, Chiapello et al. (2002) have confirmed 
the influence of the NAO on the year-to-year variability of dust export from Africa, 
using TOMS and Meteosat satellite data. 


Figure 7 shows the year-to-year variability of the NAO winter Index values from 
http://www.cgd.ucar.edurjhurrell/nao.html, and the simulated dust concentration 
at Barbados in Winter and Summer. The correlation between the NAO Index and 
the winter and summer surface concentration at Barbados are 0.67 and 0.2, respec- 
tively. For the other seasons, the correlation is as low as in summer. If one uses the 
corresponding season for the NAO Index, the correlation does not improve. Mainly 
because the pressure difference between Iceland and Portugal is much weaker, and 
thus its effects on the meteorology. 


Figure 8 shows the distribution of the correlation coefficient between the NAO 
winter Index and the winter dust emission, surface concentration and mass column, 
from 1981 to 1996. There is a high correlation (>0.9) with the dust emission from 
the Bodele depression-Lake Chad region which is the most active African source in 
winter (cf. Table 4). The year-to-year variability of surface concentration in winter 
seems to be correlated with the NAO over much of the North Atlantic and the 
western part of North Africa. Similar correlation exist for the dust mass column, 
although with a lesser extend. 


The NAO modulates the year-to-year variability of dust emission, but it is unclear 
if the year-to-year variability of dust distribution is due to the variability of dust 
emission or transport and\or removal processes. 
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6 Conclusions 


The dust size distribution in the atmosphere is simulated with the GOCART model 
from 1981 to 1996. The major features of the dust component of the GOCART 
model are that it is driven by assimilated meteorological fields and the major dust 
sources are associated with topographic lows. With these characteristics, the model 
is able to reproduce correctly in-situ and remote sensing measurements, on daily, 
seasonal and inter-annual scales. The budgets over the different continents and 
oceans have been compared. Globally, the maximum difference of annual emis- 
sion from 1981 to 1996 is about 23% with the lowest emission in 1996 (1950 Tg) 
and the highest in 1988 (2400 Tg). The dust removal by wet deposition contributes 
globally by only 10% but with important regional variations. Over the North Pacific 
the wet deposition is about 50% of the total loss. The dust deposition over the Indian 
Ocean (northern hemisphere) is as important, annualy, than over the North Atlantic. 
The annual mean North African emission is about 1400 Tg which corresponds to 
65% of the global emission. The second major dust source region is East Asia 
which contributes for 25%. Although North America has active dust sources in the 
South-west, it constitutes a sink for dust by about 30 Tg yr _1 . In North Africa and 
East Asia, there is a downward trend of dust emission of 6 and 2 Tg yr -1 , respec- 
tively. Dust emission shift from Sahel and Sahara following the seasonality of the 
ITCZ. In winter, when the ITCZ occupies its southernmost position, most emission 
is from Sahel. In summer the northward displacement of the ITCZ is accompanied 
with precipitation over the Sahel and the dry Sahara is by then the principal dust 
source. Both regions present different inter-annual variability but they both show a 
downward trend of dust emission. The inter-annual variability of dust distribution 
and emission is compared with the NAO Index. It is showed that the winter surface 
concentration is strongly correlated with the NAO Index over much of the North 
Atlantic and western Africa. We find a high correlation between the NAO Index 
and dust emission from the Bodele depression- Lake Chad region. But is unclear 
if the year-to-year variability of surface concentration is related to dust emission 
or to transport and\or removal processes. In order to differentiate the contribution 
from the different processes, it would be necessary to realize simulations with con- 
stant transport and removal processes. Also, longer simulations should be realized 
to support the strong relation between dust distribution in winter and the NAO. 
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Table 1 


Physical Properties of the Transported Size Bins 
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Fig. 2. Global Annual Budget of Dust from 1981 to 1996: Emission (line), dry deposition 
(dark gray boxes), and wet deposition (light gray boxes), in units of Tg yr -1 . 


14 



Izania (28.3N, 16.5W) Barbados (13.1 7N, 59.43W) Miami (25.45N, 80.11W) 
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Fig. 3. Comparison of the climatological monthly surface concentration, simulated (boxes) 
and observed (dots) at 12 sites in units of /jg m -3 , with their standard deviation. 
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Izana monthly concentration R=0.69 (119 points) 



Barbados monthly concentration R = 0.79 (198 points) 



Miami monthly concentration R=0.79 (184 points) 



Year 


Fig. 4. Monthly dust concentration from January 1981 to December 1996, simulated (gray 
line) and observed (iblack line) at Izana (upper panel), Barbados (middle panel), and Miami 
(lower panel), in units of ng m~ 3 . 
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Fig. 5. Comparison between the TOMS absorbing 
index (right panels) during the evolution of a dust 
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mmmmm 


0 0.500 0.600 0.700 0.800 0.900 1.00 

Fig. 6. Global distribution of the correlation coefficient between the TOMS aerosol index 
and an index calculated with the model results and TOMS viewing angles from January 
1981 to December 1990. 
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NAO index 



Year 

Fig. 7. Comparison between NAO Index (bold black) and surface concentration at Barbados 
in Winter (dark gray) and Summer (light gray) from 1981 to 1996 in units of /zg m -3 . 
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Correlation: Mass Column (Jan — Feb — Mar) & NAO Index 1981—1996 
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Fig. 8. Global distribution of the correlation coefficient between the NAO Index and the 
winter (Jan-Feb-Mar) mean dust emission (upper panel), surface concentration (middle 
panel), and mass colum (lowere panel), from 1981 to 1996. 













Table 2. Annual Budget over Six Continents in units of Tg yr- 
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Table 3 


Annual Budget over Six Oceans in units of Tg yr— 1 


Year 

N Atlantic 

S Atlantic 

N Pacific 

S Pacific 

N Indian 

S Indian 


Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 
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32 

17 

2 

76 

39 

20 

5 

126 

37 

12 

3 

1982 
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28 

15 

2 

74 

37 

20 

4 
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36 

11 
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17 
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11 
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39 

18 
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88 

45 

22 
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42 

13 
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36 

16 

3 

73 

41 

23 
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36 

13 
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30 

19 

2 

77 

42 

22 

5 
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38 

13 
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43 

20 

3 

85 

42 

25 
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35 

14 
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42 

21 

4 

85 

44 

27 

8 
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47 

15 

5 
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31 

19 

3 

77 

38 

21 

5 
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38 

12 

3 
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32 

20 

2 

70 

37 

24 

5 
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44 
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25 

20 

2 

82 

38 

24 
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38 

12 
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1992 
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30 

14 

2 
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20 
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1993 
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17 

2 

63 

32 

22 

4 

119 

30 

12 
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23 

18 

3 

86 

45 

24 

5 

140 

40 

14 

4 

1995 
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40 

16 

2 

75 

38 

23 

5 

119 

36 

12 

3 

1996 

112 

24 

18 

1 

78 

38 

22 

4 

111 

25 

12 

2 

Mean 

128 

33 

18 

2 

78 

39 

23 

5 

127 

37 

12 

3 

Dev 

12 

6 

2 

1 

6 

3 

2 

1 

9 

5 

1 

1 

Trend 

0 

-1 
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0 
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0 

0 

0 
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Table 4. Winter (January-March), Summer (June-August) and Annual Budget over Sahara (latitude >21.25°N) and Sahel (latitude <21.25°N) 
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Trend 


Table 5 

Optical Properties of Dust Particles Used to Calculate an Aerosol Index: Effective Radius 
(rjt) of the sub-bin k of the 4 transported bins (Bin), size fraction of the sub-bin (Fraction), 
Exctinction Efficiency at 380 nm (Qt), and Single Scattering ALbedo at 380 nm (a ;*). 


k 

r k 

/j, m 

Bin 

Fraction 

Qjfc(380nm) 

u^(380nm) 

1 

0.14 

1 

0.01 

0.732 

0.962 

2 

0.24 

1 

0.08 

0.276 

0.976 

3 

0.45 

1 

0.25 

3.975 

0.968 

4 

0.8 

1 

0.65 

2.427 

0.905 

5 

1.5 

2 

1 

2.354 

0.861 

6 

2.5 

3 

1 

2.228 

0.798 

7 

4.5 

4 

1 

2.182 

' 0.725 
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